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Analogous to the extension of laser cooling techniques from two-level to three-level atoms, Doppler 
cooling of molecules with an intermediate electronic state is considered. In particular, we use a rate- 
equation approach to simulate cooling of SiO"^ , in which population buildup in the intermediate state 
is prevented by its short lifetime. We determine that Doppler cooling of SiO"*" can be accomplished 
without optically repumping from the intermediate state, at the cost of causing undesirable parity 
flips and rotational diffusion. Since the necessary repumping would require a large number of 
continuous-wave lasers, optical pulse shaping of a femtosecond laser is proposed as an attractive 
alternative. Other candidate three-electron-level molecules are also discussed. 



I. INTRODUCTION 

Cold molecule spectroscopy is a promising approach 
in searches for time variation of fundamental con- 
stants [T}]3], parity violation in molecules [H [5], and 
time-reversal symmetry violation [7j. Furthermore, cold 
polar molecules have been proposed for quantum com- 
puting [21 [3] and for long-range interaction studies of 
polar degenerate gases [10] . A revievsr on recent de- 
velopments in cold molecule technology highlights these 
points, along with the difficulties associated with produc- 
ing cold molecules IT . 

Translational cooling schemes have received significant 
attention [SI [T^HTB] , and the DeMille group has recently 
demonstrated transverse Doppler cooling of SrF [TSJ [T7] . 
To date, proposed schemes have been restricted to two- 
electronic-state molecules with highly diagonal Franck- 
Condon factors (FCFs), which are molecular equivalents 
of the simplest Doppler-cooling case: two-level atoms. 
(The FCF matrix describes the coupling of different vi- 
brational levels in an electronic transition, with diago- 
nal FCFs allowing many absorption-emission cycles with- 
out vibrational repumping.) A number of two-electronic- 
state molecules have been proposed [131 [MKI^; however, 
many potential candidates have an excited state whose 
lowest level lies above the dissociation limit of the ground 
state [H]. Molecules having this characteristic suffer 
from predissociation [TH] , in which coupling between the 
excited state and an unbound level of the ground state 
leads to dissociation on time scales which can be prob- 
lematic for Doppler cooling [18 . 

Interestingly, we have found a number of molec- 
ular ions with highly diagonal FCFs, which would 
not suffer from predissociation; however, all identified 
species have an off-diagonal intermediate electronic state. 
Here, we build on existing work with two-electronic- 
state molecules and propose Doppler cooling of three- 
electronic-state molecules, potentially introducing a new 
class of coolable molecules. This development somewhat 
parallels the situation for laser-coolable atoms, in which 
repumping out of a low-lying D state considerably in- 
creased the list of coolable atomic candidates. 

One possible solution in dealing with the intervening 



molecular electronic state is to increase the number of re- 
pump lasers, in order to return population to the cycling 
transition. However, this approach would require a large 
number of repump wavelengths, because the molecule 
can decay into multiple intervening-state vibrational lev- 
els with significant probability (since the FCFs between 
the excited and intervening state are generally not diago- 
nal). Instead, we consider molecules which decay rapidly 
back into the ground state from the intermediate state. 
If the decay from the intermediate state is sufficiently 
fast, then a repump is not necessarily required, particu- 
larly for the case of ion-trap experiments in which long 
confinement times allow one to wait for spontaneous de- 
cay from a dark state. Unfortunately, spontaneous decay 
from the intervening state into the ground state results 
in diffusion of both parity and total angular momentum; 
however, the additional states can be repumped by using 
a pulse-shaped femtosecond laser, increasing the number 
of addressable transitions without a significant increase 
in the total number of lasers. 



II. GENERAL CONSIDERATIONS 

We require some number of repump lasers in order to 
scatter the 10^-10^ photons required to reach the Doppler 
limit. The total number of lasers needed to pump the 
molecule back into the cycling transition depends on two 
main factors: the FCFs, which determine the magnitude 
of vibrational branching [T31 [13 HH] , and the vibrational 
relaxation rate, which determines the magnitude of par- 
ity and total angular momentum diffusion [18j. 

As an example, we consider SiO'*' (see Fig. [T|, whose 
X^E+ and B^S]+ states have highly diagonal FCFs (see 
Fig. [2]). SiO^ also has an intervening A^H state, which 
lies 2000 cm~^ (60 THz) above the ground state. One 
would naively expect that depopulating the A state 
would require additional repump lasers. However, SiO"*" 
is representative of a class of molecules for which sponta- 
neous decays from the intermediate state to the ground 
state are on the same time scale as the intermediate state 
is populated (see Fig.[3|. 

Since the A state in SiO+ lies only 2000 cmT^ (60 THz) 



2 




FIG. 1: (Color online) Potential energy curves generated 
from the molecular constants of Cai and Francois ^20;. Life- 
times correspond to the net decay out of the v = level of 
the excited state. 



above the X state, one might expect that the B — > X 
and B — > A transitions would have similar spontaneous 
decay rates, since the transition rate is proportional to 
the cube of the transition frequency. However, the tran- 
sition moment for B — > A is surprisingly small [20] . 
while the transition moment for A — > X transitions is 
relatively large [50]. The net result is that population 
buildup in the A state is very slow. 

The low-lying A state also provides a mechanism which 
increases the vibrational relaxation rate within the X 
state. This point becomes apparent from considering 



the v' = level of the A state, which lies just below 
the v" = 2 level of the X state (see Fig. |4l. We de- 



note states belonging to the A state with 



\9v''>2,.I") 



vides a faster path to the 



The 



decay channel pro- 
states than does direct 
vibrational decay within the X manifold. This rapid de- 
cay channel results in parity flips, requiring additional 
cycling and repumping lasers, as shown in Fig. [5] Ad- 
ditionally, from the selection rule A J = 0,±1 it is clear 
that B — > A — > X decay results in rotational diffusion, 
with population building up in the \g^„ j-^^,/^) states. 

If the A state were not present, SiO'*' would behave 
similarly to SrF requiring repumping of a few vibrational 
levels of only one parity, and no rotational repumping. 
However, the addition of the A state effectively reduces 
the X state vibrational lifetimes, resulting in three sig- 
nificant changes to the repumping and cycling require- 
ments. First, the even parity states \gQ and \gQ^/.^ 
gain appreciable population and need to be repumped. 
The Lq;7i z^) and 15,^3/2) states can be driven by the 

are split 



1,1/2/ "■'^'^^ 1^0,3/2/ 

same cw source, but the \gQif.^ and \gQ^/^ 
by 129 GHz, too large a shift for typical acousto-optic or 
electro-optic modulators. Thus, the \g^ii-^ — > \^qi/2I 
and 15(^3/2) — ^ 1^01/2) must be driven by separate 
cw sources. Second, each repumped vibrational mani- 
fold {v" > 0) requires a total of four repumping wave- 
lengths (see Fig. [s]), of which only two can be derived 
from the same laser source. Third, rotational repump- 
ing of the \g^,, j//>3/2) states is required. This rota- 
tional repumping can be accomplished by driving the 




FIG. 2: (Color online) Simplified energy level diagram for 
SiO^, omitting rotational substructure and the A'^II state. 
Electronic decay from the excited state is shown, as is vibra- 
tional decay within the ground electronic state (energy not to 
scale) . 




FIG. 3: (Color online) Simplified energy-level diagram for 
SiO^, omitting the rotational substructure for all electronic 
states and the vibrational substructure of the X state. Elec- 
tronic decay into each vibrational level of the A state is shown 
along with the net decay out of the level (energy not to scale) . 
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FIG. 4: (Color online) Simplified energy level diagram show- 
ing decay to the \gQ jn) states through the intermediate j,) 
states (energy not to scale). Vibrational decay within the X 



state and decay from the \g 



p" 

v>2,.J' 



■ ) states are not drawn. 



transition (P-branch pump- 
- ^' however, each 



ing) and |<?^-->o,j">3/2) ^ 
pumped rotational level requires a different cw source, 
since the rotational splitting is 0.7 cm~^ (21 GHz). 

The actual number of driven transitions is determined 
in Sec. |III[ However, it is clear that the total number of 
required cw lasers quickly becomes unmanageable. 

An attractive alternative to using cw lasers for re- 
pumping is to use a broadband source, such as a fem- 
tosecond laser. Femtosecond-pulsed vibrational repump- 
ing has been successfully demonstrated by Viteau et ai, 
with a pulse-shaping resolution of 1 cm~^ (30 GHz) [3T]. 
With the current resolution, this technique could be ap- 
plied to rotational repuniping of the two-electronic-level 
hydrides [TS]; applying this technique to SiO'*' would re- 
quire a resolution of better than 0.7 cm~^ (21 GHz), 
which should be achievable. A typical SiO"'' repumped 
transition (e.g., v" — 1) has a lifetime of tqi — 11 fis, so 
the corresponding saturation intensity is 32 /iW/cm . If 
the frequency-doubled output of a commercially available 
femtosecond laser (P^, ~ 500 mW), evenly distributed 
across its bandwidth (A^ ~ 100 cm~^), is used for re- 




FIG. 5: (Color online) Simplified energy level diagram show- 
ing the cooling and repumping transitions for SiO"'' with ro- 
tational substructure and both parities included (energy not 
to scale). Transitions with different colors correspond to the 
addition of different lasers to the simulation, shown in Fig. [6] 



pumping, a 200 /im beam waist is sufficient to saturate 
the transition. In fact, we can take advantage of the 
comb nature of the femtosecond laser and use a con- 
trolled sweep, allowing us to periodically saturate the 
repumping transition with a larger beam waist. We note 
that for SiO^, the P branch is well separated from the 
R branch [22], reducing the complexity of the required 
spectral filtering. One could imagine driving the cooling 
transitions with cw lasers and the repump transitions 
with one or two pulse-shaped lasers. 



III. RATE-EQUATION SIMULATION FOR SiO+ 

The population dynamics of cooling are determined 
using a rate-equation approach, similar to that used 
to model population dynamics for internal cooling 
schemes [531121]. This approach allows us to demonstrate 
that population does not build in the intermediate state 
and also allows us to find the number of Doppler-cooling 
photons scattered by the cycling transition. 

Specifically, we solve 



dP 



MP, 



(1) 



where P is a vector of N energy-ordered elements, each 
corresponding to a different rovibrational state, and M is 
an N X N coupling matrix whose elements consist of Ein- 
stein A and B coefficients connecting the different rovi- 
brational states. Explicitly, the i*^ element is: 
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FIG. 6: (Color online) Total number of photons scattered 
for different laser configurations. The shaded band corre- 
sponds to cooling to the Doppler limit from room temper- 
ature. The upper limit corresponds a laser detuning fixed at 
half the natural line width, and the lower limit corresponds to 
a variable detuning that follows the temperature-dependent 
Doppler width. Cooling (Cx) and repumping (Ra;) wave- 
lengths are shown in Fig. [5] From the plot, we see that we 
require a vibrational repump for the v" < 2 levels as well as 
a rotational repump for J" = ±5/2. 



dP ^"'"^ 

— — ^ AijPi + ^ ■AjiPj (2) 

j=l j=i+l 
j=i-l ]=N 

j=l 

j=i-l j=N 
2 = 1 j=i+l 

The i*^ and j*^ states are connected by Einstein coefR- 
cients denoted by Aij , Bij , and Bji which correspond to 
spontaneous emission, stimulated emission, and absorp- 
tion, respectively, and E(ujij) is the spectral energy den- 
sity at a given frequency uJij. The Einstein coefficients 
connecting different states are calculated using Level 
8.0 |25j . The energy separations were determined from 
Ref. |26| . and potential energy curves and dipole transi- 
tion moment curves were obtained from Refs. [201 127] . 
In the three-electronic-level simulation, we include the 

IS'«"<7,J"<13/2)' I*d'<7,J'<13/2)' ^^'^ \^v<2,J<ll/2) 

states for a total of 329 states. Transitions between elec- 
tronic states are included, whereas rovibrational transi- 
tions within an electronic state are neglected, since these 
events are rare on cooling time scales. However, the effec- 
tive vibrational-relaxation rate for the X state is reduced 
since the intervening A state provides a fast-decay mech- 
anism through the X — > A i — > X decay channel, which 
is included in the simulation. 



Population is initially set to the {g^ ^^^) state, and 
the simulation is repeated for different laser configura- 
tions. In Fig. [6j we show the number of scattered pho- 
tons from the cycling transition. In Case (i), we apply 
four wavelengths to drive the \g^i^2) ^ — ^ 1^0^1/2) ^^'^ 
130^3/9) ^ — ^ 1^0^1/2) cycling transitions (lines C1-C4 in 
Fig. [5]). Parity flips resulting from the intervening A 
state require cycling both parities. Population builds up 
in the {9^1/2) ^-'■'^'^ If ?i/2) states. 

In Case (ii), we apply an additional four wavelengths 
to drive the Isj^i/j) < — > \<^Zi/2) ^"^^ \9i,3/2) ^ — ^ \^Ii/2) 
repumping transitions (lines R1-R4 in Fig. [s]). If the A 
state were not present, population would primarily go 
dark in the v" = 2 level since the v" = 1 level is re- 
pumped, and because decay into the w" = 2 level occurs 
at a faster rate than does decay into the v" — 3 level (see 
Fig. [2]). However, decay into multiple vibrational levels 
of the A state occurs on the same timescale as does de- 
cay into the v" = 2 level of the X state (see Figs. [2] and 
[3]). Subsequent decay from the A to X state is fast, and 
the FCFs are such that population primarily returns to 
the v" = level, resulting in diffusion of parity and total 
angular momentum. The resultant total photon count 
of « 3 X 10'* photons is barely above the minimum 
required number of scattered photons. 

In Case (iii), we apply an additional four wavelengths 
to drive the 15^1/2) < — > \eZi/2) and Is^g/a) < — > leji/a) 
repumping transitions (lines R5-R8 in Fig. [s]) , resulting 
in a slight increase in the photon count {N « 7 x 10^). 

In Case (iv), we apply an additional two wavelengths 
to drive the Iff^j/j) ^ — ^ |6o'3/2) branch) transitions 
(lines PRi and PR2 in Fig.[5|. Now, population builds up 
in the 19^1^/2"^ states, and the total number of scattering 
events is well above the minimum required for cooling to 
the Doppler limit (A^ « 6 x 10^). 



IV. DISCUSSION 

A. Femtosecond-assisted Doppler cooling 
candidates 

The approach of repumping by femtosecond laser pulse 
shaping increases the list of coolable molecules to include 
two-electronic-state hydrides (TH], and certain molecules 
with a third intermediate state. For example, a low- 
lying A state also occurs in molecular ions: CO"'' |28) . 
GeO+ ^29,, SnO+ ^301131], andPbO+ [31]; and for neutral 
molecules: BO [32], AlO [33], GaO [31], InO [3S], TIO, 
and N J [361 El] ■ An intermediate state that has a differ- 
ent spin multiplicity (an a state) that lies approximately 
equidistant to the excited and ground states also occurs 
in molecular ions: CH+ [38], SiH+ [39], and CC1+ [40]. 
Two cases discussed further below are CO"*", whose in- 
termediate state lies between the X and B states, and 
N J , whose ungerade intermediate state is not dipole con- 
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nected to the ungerade excited state. 

In CO"'', the B — s^X transition has a rotationless ra- 
diative Ufetime of 57.1 ns (for v = 0) with a tran- 
sition frequency of 46,520 cm^^ (215 nm), and the 
A — transition has a rotationless radiative hfetime 
of 4.03 fis (for V — 0) with a transition frequency of 
20,986 cm-i (477 nm) f28|. The radiative hfetime for 
the B — >A transition is 3.07 /xs, which we estimate by 
scahng the A — ^X Ufetime by the transition frequency 
(25,534 cm~^ (392 nm)), and by the transition moment 
which is 65% smaUer in the B — >A transition than the 
A — ^X, at the equihbrium bond length. The short life- 
time of the B — )>X transition, and the comparably long 
lifetimes of the B — ^A transition and the A — >X transi- 
tion make CO"*" a potential cooling candidate, although 
a more careful analysis of transition rates would be re- 
quired to determine the necessary repump transitions. 

Another interesting potential candidate is N^, which 
has been loaded into an ion trap via state-selective 
photoionization [JT]. Similar to SiO^ and CO"'', the 
B — >X transition of has a short radiative lifetime 
(59.1 ns for v = 0) with a transition frequency of 
25,823 cm^i (387 nm) [37], and the A — > X transi- 
tion has a radiative lifetime of 14 ms with a transition 
frequency of 9150 cm^^ (1.1 /im) [35]. However, inver- 
sion symmetry prevents decay from the B state (^S+) 
to the A state (^n„), so population does not apprecia- 
bly build up in the A state. Also, since the molecule is 
homonuclear, there is no permanent dipole moment, so 
rovibrational transitions are dipole-forbidden. Unfortu- 
nately, the FCFs of are not as diagonal as they are 
for SiO"'', so more vibrational repumping lasers would be 
required. 



B. Doppler cooling of SiO+ 

By modeling the population dynamics using a rate- 
equation approach, we determine that SiC*" is a promis- 
ing candidate for Doppler cooling. The FCFs between 
the ground and excited states are highly diagonal, with 
a B state lifetime of 70 ns. Decay from the intermediate 
A state occurs on a time scale faster than decay into it, 
so the effect of the intervening state is negligible, except 
that it leads to parity and rotational diffusion. Further- 
more, the B state minimum lies below the dissociation 
limit of the X and A states, so predissociation into the 
X or A continua [T5] does not occur. 

We must also examine the higher lying electronic states 
of SiO^. In our model, we assumed a three-electronic- 
state system. In fact, there are higher lying electronic 
states which we may accidentally excite from the B state. 
Coupling to the higher states would be of concern if they 
are unbound or if they predissociate [TS]. We are cur- 
rently in the process of performing spectroscopy with 
SiO^ to determine if these electronic states are con- 
nected via two-photon transitions from X B and sub- 
sequently from the B state to a higher-lying state. 



A' 



three-electronic-state: Cl-4, Rl-8, PRl-2 
two-electronic-state: Cl-2, Rl-2, R5-6 
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FIG. 7: (Color online) Total number of photons scattered, 
treating SiO^ alternatively as a two-electronic-state molecule 
(ignoring the A state) and as a three-electronic-state molecule 
(full treatment). The shaded band corresponds to cooling to 
the Doppler limit from room temperature for different laser 
detunings (similar to Fig. [6|. Cooling (Ca;) and repumping 
(Ra;) wavelengths are shown in Fig. [5] 



From the simulation for the three-electronic-state sys- 
tem, we find a total of four cycling transitions and ten 
repumping transitions are required. The four cycling 
transitions require three separate lasers since the 



and 



^0,3/2. 



states belong to the same spin-rotation-split 



pair (Awsr ~ 27r x 625 MHz). Of the ten repumping 
transitions, eight can be driven by the same femtosec- 
ond laser source: the Ifff 7//^i/2 3/2) — ^ 1^0^1/2)' ^^"^ 
the \3^j„^ 

=1/2,3/2) — ^ l^i^J') transitions, which span a 
region of 41 cm^^ (1.2 THz). The remaining two P- 
branch pumping transitions (15^5/2) — ^ le^f^/n)) need 



"0,5/2 



to be driven with a separate cw laser for each, or with 
another femtosecond laser. 

We also include a comparison of the number of scat- 
tered photons for SiO"*" for Case (iv) of the three- 
electronic-state system, and a simulation where the in- 
tervening state was neglected (see Fig. [7|. Here, we see 
that the intervening A state changes the nature of the 
repumping scheme; however, the A state itself need not 
be repumped. At t w 10~^ s we notice a reduction in 
counts for the three-electronic-state system, which is due 
to population building up in the |3]'^5/2) states. These 
states can be repumped using the same femtosecond laser 
present in this scheme, since the corresponding repump 
transitions fall within the laser bandwidth. 

Using a femtosecond laser does require additional con- 
sideration. Although the overall envelope of the laser is 
broad, the individual frequency components are comb- 
like. For a typical femtosecond laser, the frequency com- 
ponents are separated by 80 MHz, so it is likely that a 
repumping transition will lie within the space between 
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two frequency-comb lines. However, various modulation 
techniques, either internal or external to the laser, can 
be used to span the appropriate repumping transitions. 

V. CONCLUSION 

Similarly to how a low-lying D state does not pre- 
clude Doppler cooling of three-level atoms, an interme- 
diate electronic state does not preclude Doppler cooling 
of three-electronic-level molecules. Using a rate-equation 
approach to model the population dynamics of SiO"'' we 
determine that this species is a candidate for Doppler 
cooling of three-electronic-level molecules. The X and B 
states have highly diagonal FCFs, and population does 
not remain trapped within the low-lying A state due 
to A — !• X decay. As compared with a polar two- 
electronic-level molecule of similar reduced mass, the to- 
tal number of repumped transitions is higher as a result 
of diffusion into opposite parity and higher-lying rota- 
tional states; however, this complication does not result 
in a significant increase in the total number of lasers when 
a pulse-shaped femtosecond laser is used for repumping. 
Alternatively, if cryogenic precooling is used for trapped 
ions [3] or for neutral beams [3 [13 [IT] , the number of 
photons scattered to reach the Doppler limit is reduced 
by up to nearly two orders of magnitude. In this case, 
Doppler cooling of two-electronic-level hydrides [TS] or 
three-electronic-level oxides (see Fig. |6| requires fewer 
repumping lasers. 

For the case of molecular ions, sympathetic cooling 
by a cotrapped atomic species provides an alternate and 
quite general approach [23, 2M SI] for achieving mil- 
liKelvin molecular ion temperatures. However, direct 
Doppler cooling of the molecular ion species removes the 
need for a sympathetic coolant ion, which would be de- 
sirable in cases in which the light driving the atomic- 



cycling transition changes the molecular internal state 
or results in photodissociation. Elimination of the sym- 
pathetic coolant ions might also be desirable for cases 
in which the presence of a second species could create 
unwanted complications (e.g., quantum information pro- 
cessing using trapped polar molecular ions coupled to 
external circuits J5|). Additionally, achieving the degree 
of photon cycling necessary for Doppler cooling would 
also allow single-molecule fluorescence imaging and non- 
destructive direct detection of trapped molecular species. 

In conclusion, an intermediate electronic state is not 
prohibitive to Doppler cooling and does not require opti- 
cal pumping when population rapidly decays out of that 
state. We note that the requirements for repumping rates 
and intermediate-state relaxation rates are considerably 
relaxed in ion traps, since very long state-independent 
trapping times are routinely achieved. We have also sug- 
gested a number of other possible candidates which are 
similar to SiO"*", as well as molecules in which the inter- 
mediate state is not dipole-connected because of spin or 
gerade/ungerade symmetry. 

Currently, we are in the process of performing spec- 
troscopy of SiO^ in order to make a first observation 
of electronic states above the B state. Once the energy 
for these transitions is confirmed, we plan on proceeding 
with an attempt to Doppler cool SiO^. 
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